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vo 

O : ABSTRACT 

o 

<N 

^■T Relativistic blast waves can be described by a mechanical model. In this 

model, the "blast" — the compressed gas between the forward and reverse shocks 
C\J \ - is viewed as one hot body. Equations governing its dynamics are derived from 



conservation of mass, energy, and momentum. Simple analytical solutions are 
obtained in the two limiting cases of ultra-relativistic and non-relativistic reverse 
■^j- ' shock. Equations are derived for the general explosion problem. 

Subject headings: gamma rays: bursts — hydrodynamics — relativity — shock waves 

O 

23, 

^p H , 1. Introduction 

o 

Relativistic blast waves are believed to produce afterglow emission of gamma-ray bursts 
(GRBs). The blast wave is driven by a shell with Lorentz factor T C j ~ 10 2 — 10 3 which is 
ejected by a central trigger of the explosion. The ejecta expands and drives a forward shock 
(FS) into the external medium, and a reverse shock (RS) propagates inside the ejecta. This 
standard explosion picture has four regions: 1 — external medium, 2 — shocked external 
medium, 3 — shocked ejecta, and 4 — unshocked ejecta (e.g. Piran 2004). 

The full hydrodynamical simulations of ultra-relativistic blast waves are expensive, while 
GRB modeling requires exploration of a broad range of models. Therefore, approximate 
hydrodynamical calculations are commonly used. A customary approximation divided the 
explosion into two stages: before and after the RS crosses the ejecta. At first stage, pressure 
balance across the blast wave was assumed, i.e. the pressures at FS and RS were equated, 
Pf = p r . At the second stage, the self-similar solution of Blandford & McKee (1976) was 
applied. This description, however, has two drawbacks: (1) the approximation pf = p r 
violates energy conservation for adiabatic blast waves (by as much as a factor of 3 in some 
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example explosions that we studied), and (2) the model has to assume a sharp trailing edge 
of the ejecta, so that p r suddenly drops to zero when RS reaches the edge. A realistic model 
should describe the transition to regime p r <pj; the transition is especially interesting if 
the ejecta has a long tail. 

In this paper, we propose a "mechanical" model which does not have these problems. 
Like the previously used approximation, we assume that regions 2 and 3 move with a com- 
mon Lorentz factor T. Thus, in essence we replace the gas between FS and RS by one body, 
which we call "blast." This is reasonable since internal motions in the blast are subsonic, 
and hydrodynamical simulations confirm that T const between the FS and RS (Kobayashi 
& Sari 2000). Instead of finding T from the condition pf = p r , we leave T as a dynamical 
variable, and find its evolution from a differential equation expressing momentum conserva- 
tion. For a given instantaneous T, pressures pf and p r are found from the jump conditions 
at the two shock fronts, and the difference (pj — p r ) governs the evolution of T. 

Our model assumes spherical symmetry. However, it remains valid if the explosion is 
driven by a jet with a small opening angle 6>j Ct as long as T > Then the jet behaves like 
a portion of spherical ejecta since the edge of the jet is causally disconnected from its axis. 

In § 2, we present a useful form of the jump conditions for a shock with arbitrary 
Lorentz factor and apply them to the FS and RS of a relativistic blast wave. The equations 
of mechanical model and analytical solutions are derived in § 3. 



Ideal gas has stress-energy tensor T a/3 = (e + p)u a u@ + g a @p and mass flux j a = pu a , 
where g a/3 is Minkowski metric, u a is 4-velocity of the gas, p is rest-mass density, e is energy 
density (including rest energy), and p is pressure; all the thermodynamic quantities are 
defined in the rest frame of the gas. A shock is described by three jump conditions that 
express continuity of j a and T af3 in the shock frame (Landau & Lifshitz 1959), 

I2P2 P2 = liPi Pi, I2P2 (e 2 + p 2 ) = Tift Pi° 2 , llPi ( e 2 + V-i) + P2 = l\P\ P\<? ■ 

Here subscripts 1 and 2 refer to preshock (cold) and postshock (hot) medium, p 1 and (3 2 are 
the gas velocities relative to the shock front, and we assumed p\ — and e\ = p\c 2 . The shock 
strength may be described by relative velocity (3\2 = (Pi— /^/(l— P1P2) or 712 = (1— P^) 1 ^ 2 - 
The postshock gas satisfies e 2 /p 2 c 2 = 712, and a convenient approximation for p 2 is 



2. 



JUMP CONDITIONS 




(1) 
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It is exact if the gas is monoenergetic, i.e., particles have equal energies in the gas frame. 
Its error for Maxwellian gas is within 5%. Using equation (1), we express all quantities in 
terms of j3i 2 (or 712), which remains as the only free parameter of the shock, 

fo = P2 = 4712P1, e 2 = 47^ 2 p x c 2 , p 2 = ^ (7 2 2 - 1) pic 2 . (2) 

These equations apply to shocks of arbitrary strength, relativistic or non-relativistic. 

The blast wave has two shocks, forward and reverse. The above equations with 712 = T 
describe FS. The RS is described by the same equations when index 1 is replaced by 4 and 
index 2 is replaced by 3. Pressures pf = p 2 and p r = p 3 are given by 

4 4 
Pf = 3 ( r2 - !) Pic 2 ; = 3 (743 - 1) P4C 2 . (3) 

The approximation of common Lorentz factor of regions 2 and 3 implies that T and 743 are 
not independent. We assume that the Lorentz factor of unshocked ejecta, T e j, is known in 
the lab frame. Then one finds its Lorentz factor relative to the blast, 743 = IT C j(l — /3/3 e j)- 
For relativistic blast waves (r e j 3> 1 and r 3> 1) this relation simplifies, 

1 /r ■ r\ r 2 - r 2 

^-(f + f-)' "«=rfTP' (4) 

Thus, given the parameters of regions 1 and 4, only one parameter of the blast wave is left 
free — T. In particular, denoting p e j = p 4 , we find 



Pr _ Pcj / ej 



Pf 4r 2 p 1 Vr 2 



2 x 2 



7^-1 ■ (5) 



If the pressure balance pf = p r is assumed, it immediately determines the instantaneous T, 

'1/2 

Pf = Pr- (6) 



r ej 



i + 2r e 



^ 1/2' 
Pi X 



Pcj 



The problems of this approximation have been mentioned in § 1. In particular, energy 
conservation is not satisfied and requires a different solution. 



3. MECHANICAL MODEL 

The gas in the blast wave flows radially with 4-velocity u a = (7,7/9,0,0) in spherical 
coordinates (t,r,Q,<j)). Rest-mass conservation V M (p«^) = r~ 2 d^(r 2 pu^) = gives 
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where j- t = t| + is the convective derivative (hereafter we use units c = 1). Four- 
momentum conservation V^TJ 1 = gives two independent equations (a = 0, 1) 

1 d , o s 9/3 

~^-^ t [r h'jUa) + h'jUa— + d a p = 0, (8) 

where h = e + p. Instead of V M T M = we will use the projection u a V fl T^ 1 = which yields 
r~ 2 d fl (r 2 hu ,i ) = •ydp/dt. So, conservation of 4-momentum is expressed by 

1 d , 2 2 x <9p 2 9/3 1 d , 2 v dp 9/3 

^ ( r hl ft = ~d? ~ hl %' ^ ( r ^ = 7 dt - ^ 7 W (9) 

We apply equations (7) and (9) to the gas between FS and RS and make the approximation 

7 (t, r) = r(t), 9/3/dr = 0, r r < r < 77, (10) 

where r r (t) and rf(t) are the instantaneous radii of RS and FS, respectively. Then the 
integration of equations (7) and (9) over r between RS and FS (at t = const) yields 

^(r 2 ^) = Pr(P-Prr + \p f , (11) 

^(rW) = h r (p-[3r)T 2 +p r , (12) 
^(^r) = T^ + (h r -p r )((3-(3 r r + 3 -p f , (13) 

where £ = f* f pdr, H = f^ f hdr, P = f^ f pdr, (5 r = dr r /dt, and /3/ = drf/dt. In the 
derivation of equations (11)- (13) we made use of T 3> 1: (1) the relativistic blast is a very 
thin shell, 77 — r r ~ r/r 2 <C r, so we used 77 ~ 7> ~ r when calculating the integrals, (2) we 
used the jump condition at the FS /3/ — /3 = l/4r 2 and hf = Apf ^> pj, (3) the convective 
derivative d/dt has been replaced by (3d/ dr d/dr and T 2 /3 by T 2 in equation (12). Besides, 
the last term pj/A in equation (11) may be neglected. 

RS may not be relativistic. We derive from the jump conditions for a RS with 743 <C T, 

r ej - r 2 

P~Pr= 2r2( p2 +2r2) ' ( 14 ) 

P, = 2(^ + i-)p CJ , = p CJ , ^ = J(^ + ^ + l)p ej . (15) 

This leaves four unknowns in equations (1 1)- (13) : S, H, P, and T. One more equation is 
required to close the set of equations. We propose the following approximate relation, 

H - E = 4P. (16) 

As shown below, it is accurate in both limits of ultra- relativistic and non-relativistic RS, and 
should be a reasonable approximation in an intermediate case. 
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3.1. Relativistic Reverse Shock 

In the case of relativistic RS, 743 1, the gas is relativistically hot in both regions 2 
and 3, and obeys the equation of state h = Ap everywhere in the blast. Then X may be 
neglected in equation (16) and the condition H = AP closes the set of mechanical equations. 
Besides, we have (3 — (3 r = (2F 2 )" 1 (since r cj 3> T, see eq. 14). Equations (11)-(13) then read 

T d , 2 ^-n\ 1 1 Id, 2rr „ 2 \ o Yd, ott-^\ T 2 dH 3 3 

From the last two equations we find H(r), then substitute to the second equation and get 
the final differential equation for Y, 

dH AH A r r 

" + — = A Pl , H(r) = - / Pl r A dr. (17) 



dr r ' r 4 J 



r dr 8 \j , L J V r 4 

We have used here the jump condition p r — | (r e j/r) 2 p C j for the relativistic RS (see eq. 15). 
For external medium with a power-law density profile p 1 oc r~ k , one can simplify f p x r A dr = 
r 5 pi/(5 — k). For example, if pi oc r~ 2 (wind-type external medium), r c j = const and 
p C j oc r~ 2 , then the solution is Y = const and we find 

Pf = 3Pr, r =(^) V4 ' k = 2 - ( 19 ) 



3.2. Non-relativistic Reverse Shock 

In the case of non-relativistic RS, ^43 < 1, we have 

Pr = 4pej, Pr = \ ^Pej , At3 = 1 ~ P ~ Pr = 7^, 

6 1 cj cj 

and the equation of state u r = \p r and h r = p r + u r + p r ~ p r . The effective inertial mass 
of the blast H is dominated by region 3 as long as p r /pf 3> /3 43 . Indeed, the thicknesses of 
regions 2 and 3 are ~ r/Y 2 and ~ p 43 r/Y 2 , respectively, and the contributions to H from 
regions 2 and 3 are H 2 ~ hfjY 2 and if 3 ~ Pr^r/r 2 . This gives H 2 /H 3 ~ (p//p r )At3 <sC 1- 
On the other hand, H — H and integrated pressure P of the blast is dominated by region 2. 
In particular, P2/P3 ~ (Pz/Pr)/^ 1 3> 1- Therefore the relation H — Y, = AP remains correct 
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for blast waves with non-relativistic RS. Using this relation and subtracting equation (11) 
from (13), we derive the equation for P in the leading order of ^43, 

4 Vw«) = ?+ft- (20) 



r 2 r cj dr ' dr 



Parameters of the RS do not enter this equation. The ejecta plays the role of a piston that 
sweeps up the external medium with T pa r e j, and equation (20) in essence describes the 
shocked external medium ahead of the piston. In the case of r cj = const, equation (20) has 
the solution, 

P(r) = - r~ 8 / 3 f r 8 / 3 Pl (r) dr. (21) 
3 Jo 

Subtracting equation (13) from (12), we find HF§ = -T 2 ^ +p r (f3 - (3 r )T 2 +p r - \p f , 
which gives, in the leading order of /3 43 , 



*^=-r- (22) 



dfUs = _8 P 4 _ Pr 
dr 3 r 3 

Finally, equation (11) in the leading order becomes 

- 2Tr (rZ)=- Pe}ff (23) 

From equations (22) and (23) we find 

ps^-i^i"- (24) 

Analytical solutions can be found for ^43 and £ when the external medium has a power- 
law density profile, p\ oc r~ k , T C j = const, and p cj oc r~ 2 . Then P(r) = p x rj (11 — 3/c), 

32 PiPcjr 2 2 3 (4 - fc)rg jPl VJ _ 2 / ll-3fc \ 

' 3 (ll-3fc)(4-fc)rV P43 2 (11 - 3fc)pej' ^ 3 V 4 - A: ^ ' 1 ' 



4. DISCUSSION 

The RS may change from non-relativistic to relativistic as the blast wave propagates. 
The mechanical model formulated in equations (11)-(16) is applicable to such a general 
case. The solution can be found numerically for any given external medium pi(r) and ejecta 
T e j(t,r) and p C j(t,r). Explosions in media with inhomogeneities can be studied, as well as 
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"refreshed" explosions with an energetic tail of the ejecta. Examples of such blast waves and 
their afterglow emission will be presented elsewhere (Uhm & Beloborodov, in preparation). 

The RS shock becomes unimportant for the blast- wave dynamics when p r <C P/. This 
transition is consistently described by the mechanical model. The time of the transition 
depends on the structure of ejecta. Recent Swift observations of the plateau in the X-ray 
light curve at t Q bs < Vatcau ~ 10 3 — 10 4 s suggest that an energetic tail of the ejecta keeps 
pushing the blast as long as 10 4 s in observer's time t Q b s ~ r/2cT 2 (Nousek et al. 2006). 

The solution at late stages with p r <C Pf is obtained from dynamical equations (12)-(13) 
by setting p r = 0. In particular, for external medium with p\ oc r~ k , we find 

„, (5-k)E . . 

r =W^- p - <<p " (26) 

where E is the isotropic energy of explosion. This is close to the exact self-similar solution of 
the full hydrodynamical equations (Blandford & McKee 1976) except the numerical factor: 
T of the mechanical model is lower than the exact Lorentz factor behind FS by a factor 
of [(5 — fc)/(17 — Ak)] 1 / 2 . This difference is caused by the steep profile of the gas Lorentz 
factor 7(t,r) behind the FS at late stages. Then the approximation of mechanical model 
(d^/dr = between RS and FS) becomes less accurate, however, still gives a reasonable 
solution. At later stages, the gas motion becomes non-radial and the blast-wave equations 
must include the ^-dependence (e.g. Kumar & Granot 2003). 

The mechanical model admits a similar formulation for the case of strongly magnetized 
ejecta (Zhang & Kobayashi 2005 and refs. therein) if the jump conditions are changed 
accordingly. It can be extended to include the neutron component of the ejecta (Derishev, 
Kocharovsky, & Kocharovsky 1999; Beloborodov 2003 and refs. therein) and the effects 
of e ± -loading and preacceleration of external medium by the prompt 7-rays (Thompson & 
Madau 2000; Beloborodov 2002). 
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